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The primary purpose of the Exhauster System of the Transonic 
Turbine Test Rig of the United States Naval Postgraduate School, 
Monterey, California, is to increase the range of turbine operating 
pressure ratios by maintaining a partial vacuum inside a hood that 
encloses the turbine. The first part of this study is concerned with 
the design and the tests of the exhauster system. The experiments 
were used to establish a methpdfof predicting the operating character- 
istics of similar ejectors which is based on the universal ejection 
properties of turbulent jets. 

The second part of this study describes turbine tests utilizing 
the exhauster system that were carried out to investigate the effects 
of high values of the isentropic head coefficient and Reynolds number 
on the turbine performance. No correlation could be established be- 
tween turbine Reynolds number and performance in the range tested. 
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1. Introduction 



The Transonic Turbine Test Rig installed at the Propulsion Labora- 
tory of the Department of Aeronautics, U.S. Naval Postgraduate School, 
Monterey, California, receives its air supply from an All is -Chalmers 
VA 312 compressor. The maximum turbine inlet total pressure attainable 
from the compressor is about 43 psia. In order to increase the range 
of possible operating pressure ratios with the existing air supply, 
the turbine is enclosed in a hood which is partially evacuated by the 
exhauster system. Since the turbine can only utilize a small portion 
of the mass flow from the air supply, the remainder can be used to 
drive the exhauster. 

Very little design information is available on supersonic ejectors. 
The first part of this thesis is concerned with the evaluation of the 
design and performance of the exhauster system, and an investigation 
of the phenomena associated with a bounded, supersonic jet. 

Funds for designing and building the Transonic Turbine Test Rig 
(TTR) were furnished by the Bureau of Naval Weapons (RAPP-14). Greatly 
appreciated advice and counsel was provided by Professor M.H. Vavra. 

Mr. L.T. Clark of the Department of Aeronautics furnished a great deal 
of advice and assistance in performing experiments. 

2. Installation 

Exhauster Installation 

The complete TTR and exhauster installation is covered in detail 
by Eckert. [5] To test the exhauster separately from the turbine, a 
dummy exhauster hood was constructed for this purpose and secondary 
air was taken in from the atmosphere through a standard orifice instal- 
lation. The standard orifice installation described by Eckert was 
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equipped with a bellmouth on the intake side and a diffuser at the hood 
entrance. [4] Figure 1 shows a schematic of the secondary air flow in- 
stallation. 

In addition to the static pressure taps installed in the exhauster 
and described by Eckert, Kiel probes were installed in the nozzle 
plenum chamber and secondary air inlet pipe (Figure 1). [*5] Static 
pressure tap station designations are given in Figure 1. All pressures 
were measured in inches of mercury against atmospheric pressure as a 
reference, using 96 inch manometers with .01 inch graduations, except 
for the dr iving- nozzle total pressure which was measured on a precision 
gauge graduated in increments of .05 lbs/square inch. Secondary flow 
rate differential pressures at the metering orifice were measured in 
centimeters of mercury on a micro -manometer graduated in increments of 
.01 cm. Total temperature readings were measured with standard I.C. 
thermocouples using an ice bath as a reference. 

3. Exhauster-Only Tests 

General 

Initially the nozzle total pressure was held constant while the 
gate valve (Figure 1) was used to vary secondary flow rates. Runs 
were made with the nozzle total pressure set at 20, 25, 30, and 35 
psia. Then with the gate valve (Figure 1) fully open, nozzle total 
pressure was varied from 2 psig to the maximum obtainable from the 
All is -Chalmers compressor in increments of 2 psig. Flow rate measure- 
ments from the standard orifice, nozzle total pressure and temperature, 
and static pressure readings were taken at each operating condition. 

In order to minimize the temperature difference between the primary and 
secondary flows, the Allis-Chalmers aftercooler was run at the full 
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cold setting, giving a temperature difference of about 20 degrees Fahren- 
heit between the primary and secondary air. The tests were then repeated 
with one 40 inch long section of the mixing tube removed. 

Since the nozzle cannot accommodate the full flow rate of the Allis- 
Chalmers compressor, it was necessary to throttle the compressor flow, 
discharging a portion of it to the atmosphere. The maximum total pres- 
sure obtainable at the driving- no zzle was limited by approaching temper- 
ature and surge limits on the compressor. A six foot flexible pipe 
connection between the driving- no zzle and air supply (Figure 1) pro- 
duced an unexpectedly large pressure loss (about 8 psig static pressure) 
limiting the nozzle total pressure to a maximum of 36 psia. 

The entire system was allowed to stabilize before taking readings 
at each point of operation. The system was very stable except at high 
secondary flow rates, where the static pressures at stations one through 
four showed a large low frequency variation of as much as 2 inches of 
mercury. However, this instability was not reflected in the total 
pressure maintained inside the hood. The instability in the mixing 
tube at high flow rates is contrary to both experimental evidence and 
theory given by Alexander, Baron, and Comings. [3] These authors found 
that the instability in the mixing tube occurred at zero secondary flow 
rates and hypothesized that this was due to the natural function of the 
jet in entraining air from its surroundings, requiring a reverse flow 
eddy in the absence of any secondary flow. At zero secondary flow 
rates, the TTR exhauster system was very stable, even at the maximum 
driving nozzle total pressure ratio. It is possible that hood leakage 
may have furnished a slight amount of secondary flow to stabilize the 
jet at zero secondary flow rates, although no evidence of any leakage 
was detected. 
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The first attempt to obtain velocity profiles was made with a W-157 
probe built by the United Sensor Corporation, East Hartford, Conn., which 
is a wedge shaped probe with an included angle of about 34 degrees. The 
probe has a total pressure tap located in a small cut-out section and 
static pressure taps located on the face of the wedge. A mercury man- 
ometer was used to set the probe to a zero yaw angle by equalizing the 
static pressures from the two face taps. 

Although it was realized that the detached shock Mach number for 
the wedge would be above the Mach numbers encountered in the flow, it 
was hoped that a calibration curve could be obtained over the Mach num- 
ber range of interest (.7 to 1.9). The supersonic wind tunnel with 
variable Mach numbers available at the U.S. Naval Postgraduate School 
is an AMRAD Model W-4 miniature blowdown tunnel. Since the physical 
dimensions of the probe precluded use of the actual probe for any cali- 
bration work, a full scale model of the probe tip complete with total 
and static pressure taps was constructed so that it could be mounted 
in the tunnel. Attempts to calibrate the probe on two occasions were 
completely futile since none of the readings obtained were repeatable 
within a tolerance of t 50 percent. 

Flow over the wedge was examined visually by means of a Schlieren 
system but failed to show any unusual conditions which might have ac- 
counted for the wide variations in dynamic pressure. It is possible 
that the tunnel calibration curves are inaccurate or that repeatable 
Mach numbers cannot be obtained with the tunnel. The probe was tested 
for leaks, and none were found upon completion of testing. 

While waiting for the construction of the probe model, velocity 
traverses were taken with the W-157 probe with the concept that flow 
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yaw angles could be checked even if the probe could not be calibrated. 

No evidence of any large swirl components was found since the yaw 
angles did not exceed 5 degrees and were usually less than 1 degree. 
However, upon completion of the traverses, it was found that the W-157 
probe had developed a leak between the total and one static tap; there- 
fore, no conclusions can be based upon data obtained with this probe. 

The next attempt to take velocity traverses was made with a pitot - 
static tube 1/16 inch in diameter mounted on a \ inch stainless steel 
tube. This tube produced very erratic readings attributable to severe 
probe vibration. Since the probe did indicate a very unsymmetrical flow 
pattern, this possibility was checked by reading static pressures from 
one orifice at a time at each station. All of the pressure taps pro- 
duced readings within .05 inches of mercury of each other showing no 
indication of unsymmetrical flow. 

Since the difficulties seemed to be due to probe vibration, a 
total head probe (Figure 2) was constructed so that it could be sup- 
ported from both ends. The idea of constructing a pitot-static probe 
was discarded initially, since the holes through which the probe had 
to be inserted were too small for a static probe with the orifices lo- 
cated far enough back of the tip to conform to standard practice (eight 
to ten diameters). Therefore, it was assumed that the static pressure 
was uniform across the section and total head traverses were made only. 
The total head traverse at station one indicated a lower total pressure 
at the centerline than at the off-center stations. Since this phenomena 
was not expected at the time, a small non-standard static probe was 

c 

constructed anyway (Figure 2) to take static traverses. The velocity 
profiles shown in Figures 3 to 5 were obtained with these measurements. 
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No attempt was made to calibrate the probes due to lack of facili- 
ties. In addition, the probes could not be inserted or withdrawn while 
the system was in operation without damage to the probe, so that a 
slightly different equilibrium point was attained for each station. 

Data Reduction 

Secondary flow rates were measured by the orifice installation in 
the secondary flow inlet. Calculation of flow rates was accomplished 
by computer program based on the orifice formulas used by Eckert. [4] 
Agreement of flow rates calculated from vena contracta and flange tap 
differential pressures were within 2 percent for all flow rates and 
generally within 1 percent or less. 

Nozzle flow rates were calculated from assumed polytropic expo- 
nents and measured values of nozzle total temperature and pressure. 
Calculations for nozzle flow rates in the subsonic region utilized 
isentropic flow and one-dimensional normal shock relations. 

The velocity profile data reduction was accomplished by means of 
a computer program using the Rayleigh formula, with the known local 
static and total pressure behind the shock. Standard compressible 
flow relations were used for subsonic velocities. 

4. Exhauster Performance 

General 

The classical method of analysis is based on a one- dimensional 
momentum analysis; however, if nothing is known about the secondary 
flow except its initial total pressure and temperature, the problem 
cannot be solved since the number of unknown quantities exceeds the 
number of independent flow equations that can be formulated. There- 
fore, it became necessary to find an additional relationship between 
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the primary and secondary flow rates. 

Of all the various approaches dealing with the problem of turbulent 
jets, the one that appeared to have the greatest potential is a hypo- 
thesis set forth by 0« V. Yakovlevskii. [l4] His hypothesis, based on 
a wide survey of experimental and theoretical data^ is that in the dis- 
charge of a gas jet, the ejection properties are the same as those of 
a geometrically similar submerged jet of an incompressible fluid with 
the same initial impulse. A submerged jet is defined as a free jet 
discharging into a fluid at rest. This hypothesis leads to the con- 
cept of the universality of the ejection properties of a turbulent jet 
regardless of the conditions of discharge. 

Assuming that the basic property of a jet is its impulse, and that 
the basic characteristic of such a jet is the ejection or entrainment 
of mass from the surrounding media, an arbitrary jet and a standard, 
isothermal, submerged free jet discharging from geometrically identical 
nozzles with identical impulses in the initial cross section should 
have similar ejection properties. 

The basic part of Yakaievskii' s analysis is presented here as a 
basis for further development, retaining the author's notation: 

Restating the basic hypothesis in mathematical form, 

m dG; r d Gc^ 

dx JT 

where: 

G = gas discharge, pVA in the initial cross section 
x = longitudinal coordinate in the direction of flow 
i = index referring to an arbitrary jet 

a = index referring to a standard, submerged, isothermal jet 
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Then, the relationship between the standard jet and the arbitrary 



jet is expressed by the ratio, 




o 

where is the density of the surrounding medium and 9 is the density 
ratio. Let: 



(3) 




where G 0 is the gas discharge at the initial cross section and, 



(4) 




X 



r 0 

where r Q is the initial radius of an axisymmetrical jet. 
now be written in a non-dimensional form: 



Eq. (2) can 




This hypothesis amounts to the assumption of the universality of trans- 
verse and longitudinal velocity profiles. 

It is shown by Yakovlevskii that the increase in flow rate in some 
arbitrary section of the submerged jet (Figure 6) is, [l4] 



( 6 ) 



dG * 2 tr 9 



where : 



9^= density in the undisturbed part of the flow 

Voo = instantaneous velocity in the y direction in the undisturbed 
part of the flow 

Putting Eq. (6) into a non-dimensional form, 



(7) 



O 



dG . 

cU° 



2 S. Uvv> n° Voo 
Vo ‘-'o “ Uv* 
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where : 



9 = 9*° = the ratio of the density of the undisturbed portion 

p o of the surrounding flow to the density of the jet 

T = radius at some arbitrary point of the jet 

= instantaneous centerline velocity in the x direction 

Voo = instantaneous velocity in the x direction in the undisturbed 
part of the flow 

Mo = instantaneous velocity in the x direction of the jet discharge 
He goes on to state that the following relationship is valid for this 
case: 

T Mwi _ i 

r * u ; 

and that for small values of Q , is a constant. Thus, 



(9) 



- 4= z const 

dx Vp Uw> 



Then, using Eq. (7) he gives the formula for the initial cross- 
section of an axisymmetric jet as, 



( 10 ) 



G 



° °< i (x° +61) 



where ©<, and 0, are experimentally determined constants. The follow- 
ing values are given by Yakovlevskii for the case of an incompressible 
submerged free jet: [14] 

<X| = 0.11 to 0.20 
£1 = -4.3 

The minus sign attached to o< 4 is due to the Russian practice of basing 
jet coefficients on the assumption that the jet is issuing from a point 
source located at some distance in front of the nozzle^as Eq . (10) be- 
comes meaningless near the nozzle. Taking the optimum length of the 
mixing tube as 8.5 nozzle diameters, an average of the commonly used 
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value of 7 to 10 diameters, and taking <X| as 0,115, there is for an 



4.76 = driving nozzle diameter (inches) 

8.5 = optimum length of mixing tube (driving nozzle diameters), 
taken as the value for x 

4.76/2 = initial radius of the axisymmetrical jet, r Q (inches) 



The optimum length chosen for the mixing tube was verified by re- 
moving one 40 inch section of the mixing tube. The tests showed that 
almost ‘identical operating characteristics were obtained for the 40 
and 80 inch long tubes. Comparing the coefficient of discharge of Eq. 
(11) with the experimental values of the present data (Table I) shows 
very good agreement between the two provided that a suitable choice is 
made for o(, . However, on the basis of information available, and due 
to the wide variation of o<i from author to author, the arbitrary choice 
of c*i to fit experimental data cannot be defended. It is noteworthy 
though that the term'G° is nearly a constant throughout the entire range 
of operation up to the point of secondary flow choking. Therefore, an 
attempt was made to place the flow relationship of Eq. (10) on a firmer 
basis. Using the standard definition of impulse function, 



the nozzle impulse function was plotted versus where p o is the 



density of the nozzle discharge. Figure 7 shows that is a linear 



function of the nozzle impulse function. If the impulse function is 



assumed 




of unity, 



(ID 




where: 




from about 1.2 to 0.96 for the present tests. 



(13) 



F = PA 
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plotted versus the secondary flow rate (Figure 8), a linear relation- 
ship is evident up to the point of secondary flow choking, and thus the 
secondary flow rate for this particular system at a given nozzle im- 
pulse can be determined from Figure 8 for unchoked secondary flow* The 
linear relation holds only for the case where the nozzle discharge is 
supersonic. Since the geometry of the nozzle is fixed, there was no 
way to test the generality of the coefficient of discharge determined 
from Figure 8; however, since the basic hypothesis appears to be true, 
it should provide a reasonable first approximation for a geometrically 
similar system. 

The reason that the linear relationship between the nozzle im- 
pulse function and the secondary flow rate does not hold for the sub- 
sonic case is probably due to the fact that one-dimensional isentropic 
and normal shock relations were used to calculate nozzle exit condi- 
tions, and these are not truly representative of the actual flow 
phenomena. 

The next problem was to find some means of determining the required 
mixing tube length. Since the original hypothesis was based on the as- 
sumed universality of velocity profiles, it was assumed that a determi- 
nation of the necessary length of the mixing tube could be based on 
incompressible flow. However, even with this simplification, the prob- 
lem was not tractable due to the presence of a longitudinal pressure 
gradient (Figures 9, 10) and the problems associated with turbulent 
flow. Therefore, an empirical approximation was sought. 

Viktorin carried out a series of experiments in 1941 using water 
as the driven and driving fluid. [12] His experimentation was based on 
an earlier theoretical analysis performed by Flugel. Viktorin found 
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that the velocity distribution in the mixing chamber behaved in the same 

1/3 

fashion as a circular wake with the width proportional to about x and 

-2/3 

the centerline velocity proportional to about x . The assumption is 
made that the mixing process reaches an optimum point when the core of 
the jet contacts the mixing tube, and that any further increase in mix- 
ing tube length will result in a trade off between increased mass flow 
rate and increased mixing and frictional losses. Assuming that the 
primary fluid emanates from a point source located at some point up- 
stream of the nozzle in order to find a constant of proportionality, 
it was determined that the mixing process would reach an optimum state 
about 40 inches downstream of the nozzle. This was verified by the 
performance data (Figures 12, 13 and 14) for the system with one 40 
inch section of the mixing tube removed. 

A comparison of predicted centerline velocities with experimental 
data (Table II) shows that the predicted values are within 10% of the 
experimental values for stations two through four. The discrepancy at 
station six can be attributed to the fact that at this point the flow 
is completely mixed and is now characteristic of turbulent flow in a 
pipe. It is thought that the flow at station one differs from the pre- 
dicted value due to off-design operation of the nozzle, i.e., the 
nozzle discharge static pressure is greater than the receiver static 
pressure, and the nozzle never reached the design point during the 
experimentation. Examination of the total pressure traverse shows 
that the maximum total pressure was not obtained on the axis of the 
jet. This matches experimental data given by Abramovich. [1] He shows 
that the jet actually has a small subsonic region along the centerline 
and that the centerline total pressure is less than the off-axis total 
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pressure until the flow has become sufficiently adjusted to become iso- 
baric. 

The analysis up to this point has considered only the case of un- 
choked secondary flow. As is evident from Figure 12, the secondary 
flow rate becomes almost independent of the nozzle total pressure above 
a pressure ratio of 2.3. A supersonic ejector has two regions of op- 
eration. One region, the so-called mixed region occurs when the driving- 
nozzle total pressure ratio is sufficiently low to make the ejector flow 
dependent on the back pressure. If the driving-nozzle pressure ratio is 
increased sufficiently, the flow will become independent of the back 
pressure, since it becomes supersonic. As the nozzle pressure ratio is 
increased even further, the secondary flow will eventually choke at the 
point of minimum area. 

The supersonic region is characterized by the condition that the 
secondary flow static pressure at the driving nozzle exit plane is less 
than the nozzle exit static pressure. This point of operation was not 
reached with this particular system, and the system operates in the 
mixed region. However, since secondary flow choking does occur, the 
possibility of sonic flow at the minimum secondary flow area was checked 
over a wide range of nozzle pressure ratios, and it was found that the 
secondary flow does not choke at this point. Calculation of the poly- 
tropic exponent and the flow function (Equation 14) between the secondary 
flow inlet and the nozzle annulus (Table III) shows that at the higher 
flow rates, the process is nearly isothermal and that the flow function 
approaches the critical flow function (Equation 17) indicating choking 
at some point between the inlet and the nozzle annulus. 

If the secondary flow is not restricted prior to reaching the 
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nozzle exit plane, flow characteristics can be found using a method 
developed by Fabri. [7] A more elaborate method utilizing the method 
of characteristics was presented by Addy. [2l Both of these methods are 
concerned with flow in the supersonic region. For flow in the mixed 
region, both authors use a form of the one-dimensional momentum anal- 
ysis. Hovever, as mentioned before, solution of the momentum equation 
depends on control of the secondary flow rate or pressure ratio. 

The empirical method of determining the ejector pumping character- 
istics has some serious drawbacks as presented here. The coefficient 
of ejection which is based on empirical data is an unknown function of 
system geometry, and since the viscous effects have been ignored, it 
does not provide any means of predicting the pressure distribution in 
the mixing tube. Figures 9 and 10 show the mixing tube static pressure 
profiles. No theory could be developed to produce the mixing tube 
static pressure distribution. 

Momentum Analysis 

Once the values of secondary flow rate are determined from Eq. 

(10), a one-dimensional momentum analysis can be conducted to determine 
exhauster performance. The analysis was designed to parallel the pre- 
liminary design calculations by Vavra in order to determine the validity 
of the on®— dimensional approach. [lO] In order to simplify the equations, 
polytropic process relations as well as the non-dimensional flow func- 
tion^ were used throughout. These were taken from Vavra. [9] 

For an expansion process, 

[(&.)* -(&>*] 
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where : 



( 15 ) 



9 - 

^ 1 i - y\p(^-^ 



n = polytropic exponent 
P = ratio of specific heats 
W = flow rate (Ibm/sec) 

T t = inlet stagnation temp. (°R) 

o 

A = flow area at station of interest (in^) 

P to = inlet stagnation pressure (psia) 

P = static pressure at station of interest (psia) 

R = gas constant for air (ft Ib/lbm °R) 
g' = universal gravitational constant (Ibm-f t/lb-sec ) 
For a compression process, 




where the subscripts 1 and 2 refer to the inlet and station of interest 
respectively. Figure 1 shows the exhauster with the station designa- 
tions used in the analysis. 

For the exhauster operating conditions the flow through the nozzle 
will be choked. For this condition, the flow function,<£, becomes 



( 17 ) 





It was assumed that the nozzle total pressure and temperature are 
known. Assuming a polytropic exponent for the flow in the nozzle to 
the throat, can be calculated. Then, another polytropic exponent 
was assumed from the throat to the nozzle exit so that, 

(is) 
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Expressing the dimensionless flow function in terms of pressure ratios, 



<19) *. = 

Since <J) c is known from the choked flow condition, and and T^ are 

known, it was possible to solve for P IN by iteration. Having found 



( 20 ) 


T '~ = T ™ (&) " 


( 21 ) 


V ) N - ^gTCpHTN-Tlf^ 


( 22 ) 


Wn 5 fePTNANW RfrN 



Thus the complete operating conditions have been determined for 
the driving nozzle. It was assumed that the air leaving the nozzle 
would undergo a full expansion without oblique shocks in order to con- 
tinue with the solution. If the nozzle total pressure is low enough, 
the back pressure at the exit plane will be high enough to produce 
shocks in the nozzle with attendant subsonic flow at the nozzle exit. 

An investigation conducted on the experimental data showed that the 
flow from the nozzle would be independent of the back pressure whenever 
the nozzle total pressure exceeds 22 psia. Therefore, this particular 
analysis is only valid for nozzle total pressures above 22 psia. 

Secondary air flow enters the flow measuring orifice with a total 
pressure, P ^q, and a total temperature, T.j,q, equal to atmospheric con- 
ditions. Since the process is essentially adiabatic, T^q will remain 
a constant. Then for the flow from the hood to the nozzle exit plane, 



( 23 ) 
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Because the weight rate of flow of secondary air is known for a given 
set of nozzle conditions, a ratio of Pit^TD can assumed, along with 



a value of PjT* For assumed pressure ratio, 



(24) 



T.T **Ttd (Ql) 



Oil 

n 



(25) V I T =^;>g J c p (Ttd - Tit ) 

Writing the momentum equation between stations 1 and 2, by assuming 
uniform velocity profiles, uniform static pressure distributions, and 
equal static pressures at the nozzle and secondary air discharges, 

(26) ( Wn Wr jV?t _ Wn Vjn _ Wt Vit = A2i(Phm-P 2 0«. F *f 

3 3 3 



The friction force Ff, was expressed as a pressure loss: 



(27) F-f : 1L Dl & P -A2\ a P 

H 



where, 

(28) AP - i V* 

Di 2 

L = length of duct 
D 2 = duct diameter 
V = average velocity 
p = average density 
f = friction factor 

f is a function of Reynolds number and pipe roughness. With an 
approximate range of Reynolds numbers from 1. 7-2.2 x 10^ and a surface 
roughness of .00015 for commercial steel pipe, a roughness to diameter 
ratio of .0024 was obtained. Entering a Moody diagram, a value of 0.017 
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for f was obtained. Then, a aP was calculated and subtracted from P_._ 

IN 

to obtain P*. 

Rewriting equation (26) , 

(29) ^ Wn +Wt JV?i - WnVin- WtVit - An (P*-P 2 i) 

The total temperature after mixing was taken as a weighted average: 



( 30 ) Tt 2 : WtTto Wn Ttn 

Wn + Wt 

With Wn -Wt - W and writing the continuity equation, 



(31) V21 -- WRT21 : WR /Tt 2 . - V2* \ 

A21 P21 P21A21 \ 2 Tcp / 

Solving (31) for P^ and with (30), 



(32) P 2.1 - W R / Wt Tt o ♦ Wn Ttn _ Vat ^ 

An V21 [ W icj 7 Cp I 



Combining (30) and (32), 



( 33 )/ Wn + Wt V21) -Wn Vin - Wt Vit =Aai P *- WR 
\ 9 / 9 9 V21 

2 

/ WtTto * WnTtn - V2 1 \ 

\ W igltp / 



Multiplying (33) by and expanding, 



(34) V 2 | 2 f WiM 1 Wt „ WR ) - V2I ( Wn Vin + Wt Vit 

\ 9 2 3 Tcp / ' 9 

TA 21 PM + WRT * O 



Let: 
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(37) 



B 3 : WRTn 



Solving the quadratic (34) for velocity, 



(38) 



V21 -- - 82. +\// 82^ • B 3 
2 61 ~ Vl^- Bi / Bi 



This equation yields both a subsonic and a supersonic solution; 
however, only the subsonic solution was pursued since the supersonic 
solution did not provide sufficient pressure recovery under the opera- 
ting conditions investigated. With the known quantities, 



(39) 



(40) 



T 21 1 Tt 2 - \h. 1 

2-cj Jc.p 

Pt 2 = Pz> | T T2 ^ 



fr-l 

r 



For the subsonic solution, the flow undergoes a compression pro- 
cess in the diffuser: 



<41) $ 

This equation was solved by iteration to find P3. If P3 does not match 
P AMB (atmospheric pressure) a new value for P^ N was assumed and the 
process repeated until P3 matched 
Having obtained P3, 



(42) 



T 3 = Tz \ 



(&) 



nc ~ L 
nc 



(43) 



(44) 
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Test Results 



The equations just described were used to predict exhauster per- 
formance utilizing a computer program. Since the primary purpose of 
the exhauster is to maintain a partial vacuum inside the hood, this was 
the major point of interest. A comparison of experimental and theoret- 
ical data is given in Figure 11. The discrepancy is due to the fact 
that the one-dimensional analysis ignores conditions from the inlet to 
the hood, assuming that the polytropic exponents are constant. This 
is not the case since secondary flow conditions are dependent on the 
varying polytropic process from the inlet to the hood. Table III shows 
the variation in polytropic exponents. 

For the case where the secondary flow into the hood is a constant 
(i.e., turbine operation), the one-dimensional momentum analysis ap- 
pears to provide a good approximation (Figure 11). It is to be noted 
that the experimentation was conducted with a nozzle total temperature 
about 20°F above the secondary flow total temperature. Increasing the 
nozzle total temperature will result in a very slight decrease in ex- 
hauster pumping capability, and for the allowable range of turbine 
operating temperatures, exhauster flow is virtually independent of 
total temperature. 

5. Exhauster Design Evaluation 

General 

The most recent work on the design of ejectors available was pub- 
lished by Engle in 1963. [6*] He presents a set of equations to be used 
to optimize the design geometry based on a one-dimensional analysis of 
incompressible flow. The present exhauster system meets the criteria 
set forth by Engle for driving-no zzle placement and secondary flow 
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entry. For optimum operation, the two streams should be parallel upon 
mixing tube entry to minimize losses. Secondary flow static pressure 
profiles (Figure 14) show a continuous decrease in static pressure up 
to the nozzle exit plane except at very high nozzle pressure ratios. 

The increase in static pressure at station 28 at high nozzle pressure 
ratios is probably due to separation. 

At high secondary flow rates, the mixing tube static pressures 
show a fluctuation of 1.0 to 2.0 inches of mercury. However, the flow 
rate and hood total pressure remained steady in this regieme of opera- 
tion. Since the secondary flow is choked at these pressure ratios, it 
was thought that the secondary flow rate might not be sufficient to 
stabilize the jet. This was not the case though for this system. The 
exhauster operation was very stable at high nozzle pressure ratios with 
no secondary flow. The instability only appears at high pressure ratios, 
without being felt inside the hood. The probable cause of the instabil- 
ity is the formation of a recirculation eddy due to separation effects 
along the mixing tube since the secondary flow is under an adverse 
pressure gradient. 

Engle riO] also presents a method for determining the required area 
of the mixing tube. This is based on continuity, and is, 

(45) A -- W 

Vai ViN 9 
Vim 

Substituting experimentally determined values into this equation for a 
nozzle pressure ratio of 2.22, 

(46) A • (8.011*3.3) -.686*+* 

7oo(.M4 1 ) IToTTTj 
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where £ is an average density based on the average velocity computed at 
station five. The actual area of the mixing tube is .306 ft^. Since 
Engle states that in practice the mixing tube area should be some 30 to 
50 percent larger than the calculated value, the mixing tube area of 
this system is about 30% of the desired optimum. 

Optimum mixing occurs according t<5' Engle when the ratio between 
the velocity in the core and in the annulus is between .7 and .8 and 
that the length required is a function of the nozzle to secondary flow 
area. £6} In general, the optimum mixing length is between 7 and 10 
diameters of the mixing tube. Removal of one section of the mixing 
tube left the exhauster operation virtually unchanged except at the 
higher pressure ratios. The velocity profiles taken at the diffuser 
entrance (Figure 5) for the short mixing tube indicate that the ratio 
between the core and annulus velocity increases with increasing nozzle 
pressure ratio, and falls short of the .7 to .8 desired velocity ratio 
which accounts for the reduction in pumping capability at the higher 
pressure ratios. 

Optimum mixing tube length for this particular exhauster falls 
between the short and long tubes. However, as indicated from the ex- 
hauster performance, a wide variation in mixing tube length has little 
effect on system performance, since varying the length to diameter ratio 
from 5.46 to 10.9 had only a small effect on the operating characteris- 
tics . 

It was not possible to duplicate turbine operating conditions as 
far as secondary flow rates and nozzle total pressure are concerned. 
However, extrapolating from experimental data, it should be possible 
to maintain a pressure ratio of about 5 to 1 for the converging -diverging 
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stator and a pressure ratio of about 6.3 to 1 for the converging turbine 
rotor iator for the transonic turbine. Exhauster operating conditions 
can be determined from Figures 11, 12, and 13. 

Discussion 

The basic hypothesis set forth by Yakovlevskii appears to be valid 
since the ejection property of a given nozzle is a function of its im- 
pulse, and empirical relations found for incompressible flow can be 
used to find a first approximation for the case of compressible flow. 

If the secondary flow rate or pressure is known, a one-dimensional 
momentum analysis can be used to determine operating conditions pro- 
vided that the secondary flow is not choked. If the secondary flow 
is choked at the mixing chamber inlet, the method of Fabri can be used 
to predict operating characteristics. The method based on the hypothe- 
sis of Yakovlevskii has the advantage of producing a good approximation 
of the pumping characteristics with relative ease. However, the gen- 
erality of the relationship between the secondary flow and nozzle im- 
pulse function was not tested for other ejector configurations. 

Determination of the required length of mixing tube for a given 
ejector is still largely a matter of experience \ however, it was shown 
that a wide variation in mixing length had very little effect on system 
operation. Mixing tube length to diameter ratios of from 7 to 10 
should produce optimum operation. Calculation of the required mixing 
tube diameter based on the equation given by Engle indicates that the 
mixing tube diameter falls short of the desired optimum. 

Although unstable conditions did appear in the mixing tube at high 
flow rates, the instability was not reflected in the hood total pres- 
sure. This fact coupled with the extrapolated performance figures 
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indicates that the exhauster will provide turbine operating pressure 
ratios in the design range. 
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Part II 

Tests of a Reaction Turbine 
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1. Introduction 



The second portion of the thesis describes a turbine test program 
carried out with the exhauster system to investigate the effects of 
high values of the isentropic head coefficient on turbine performance. 
The turbine tests were used also to correlate turbine performance with 
Reynolds number. 

2. Turbine Installation 

The turbine installation is identical to that described by Eckert 
except as noted here. [5] The turbine is an ARES MOD II instead of the 
Transonic Turbine. A cross section of the turbine is shown in Figure 
15. Figures 16 to 18 show details of the rotor and stator construction. 
The instrumentation was modified by placing additional pressure taps 
in the shroud as shown in Figure 15. . It is to be noted that the pres- 
sure tap locations have been projected into the one dimentional view 
while in actuality they are spaced around the periphery of the shroud. 
For this particular test program, the axial clearance was set at 0.41 
inches while the radial clearance was set at 0.033 inches. 

The cover plate flexure instrumentation mentioned by Eckert was 
installed to permit measurement of the torque and the axial force on 
the cover plate. [5] ' 

3. Turbine Test Program 

Flow Nozzle Calibration 

The turbine flow nozzle installation and calibration techniques 
used are covered by Eckert. W Early calibration runs conducted by 
Eckert and Mr. L. T. Clark indicated that at low flow rates, the nozzle 
coefficient seemed to be a function of supply total pressure. Since 
this violates the principle of similarity, additional calibration runs 
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were made at supply total pressures of 20, 22, and 24 psia. Data was 
reduced using computer program FLOCAL with the exception that the ex- 
pansion factor for the nozzle was changed to that given by the ASME 
Power Test Codes, [s] Figure 19 shows the nozzle coefficient as a func- 
tion of Reynolds number. The nozzle coefficient was based on vena 
contracta tap data since these have a lower tolerance. Examination of 
earlier calibration runs reveals a fairly large scatter (as much as ten 
percent) between vena contracta and flange tap data at low flow rates. 
For the calibration run shown in Figure 19, the difference between vena 
contracta and flange tap calculated flow rates is less than one per- 
cent. Therefore, a nozzle coefficient based on Figure 19 was used. An 
analytical expression for the flow nozzle coefficient as a function of 
Reynolds number was found by using the method of least squares to ob- 
tain a polynomial approximation. This was accomplished by a computer 
program yielding an eighth order polynomial for the best fit. The poly- 
nomial is, 

C = -69.4248 + 65.42325x - 25.4166x 2 + 5.38697x 3 

(47) 

-0.67816X 4 + 0.051422X 5 - 0.0022472x 6 + 0.000049356x 7 
-0.000000 3595 9x 8 

where x = Reynolds number multiplied by 10 -3 . The maximum deviation 
in the flow range of interest (1.6 to 4.5 Ibm/sec) occurs at a Reynolds 
number of 6 x 10^ and is .003. The polynomial is only valid for Rey- 
nolds numbers in the range of 4 to 12 times 10^. 

Turbine Flow Rate Determination 

The turbine flow rate is equal to the total rate measured by the 
nozzle less the plenum labyrinth leak rates. Total nozzle flow was 
determined by means of the ASME standard equations using an initial 
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nozzle coefficient of 1.03. [is] Vena contracts tap data was used 
throughout, with differential pressure measured on a 96 inch water 
manometer, and flow nozzle static pressure measured on a 96 inch mer- 
cury manometer. Then, after calculating the Reynolds number, the flow 
nozzle coefficient was corrected by application of the eighth order 
polynomial described above and a new flow rate determined. 

Calculation of the plenum labyrinth leakage flow rates utilized 
the expressions given by Eckert which were determined from experimental 
data. [4] . 

Turbine Tests 

The turbine was rotated on 14 November 1966 with the bottom portion 
of the hood installed to check the instrumentation. The shaft seal 
labyrinth was bored out to 1.258 inches diameter to accommodate a* new ' 
quill shaft without recutting the chambers for a seal tooth clearance 
of 0.040 inches. The turbine was run up to about 13,000 RPM. Upon 
shutting down, it was noticed that, the quill shaft had rubbed on the 
labyrinth, scoring the shaft. The shaft was buffed and reinstalled 
with the labyrinth realigned. ' When the turbine was rotated again on 
16 November, shaft rubbing occurred at. about 7,000 RPM. The shaft was 
measured and found to. be 0.005 inches out of round. It was then turned 
down 0.005 inches and the labyrinth enlarged so that after alignment 
the clearance was 0.030 inches radially at the dynamometer end and 0.015 
inches radially at the hood end. This prevented further rubbing when 
the turbine was rotated on 18 November 1966 with the full hood and ex- 
hauster installed. The present labyrinth is only a temporary device 
and will be replaced. Since the only result of the increased leakage 
into .the hood is the degradation of the system's capability to maintain 
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a partial vacuum inside the hood, the existing labyrinth is suitable 
for continued testing until a replacement is obtained. The purpose of 
the first run was to investigate the performance limits of the system. 
With the ARES MOD II turbine installed, a pressure ratio of 1.79 was 
attained. Since the turbine efficiency is much greater than that of 
the exhauster, a reduction in turbine flow rate with the corresponding 
increase in exhauster flow rate resulted in a decrease in pressure ratio. 
The exhauster is capable of maintaining a partial vacuum of one-half an 
atmosphere in the hood with a turbine flow rate of 3.88 Ibm/sec. Since 
this particular turbine has been the subject of an extensive test pro- 
gram without the hood installed reaching a maximum pressure ratio of 
1.55, it was decided to run at pressure ratios of 1*3, 1.4, 1.5, 1.55, 
and 1.65 with values of the isentropic head coefficient varying between 
the design value of 2.609and 4.5. 

Turbine testing was started on 22 November 1966. During the run, 
the RPM readout was checked by means of a strobotac and it was found 
that the RPM read out was faulty. The system was checked thoroughly, 
and the difficulty finally corrected by replacing the leads between the 
counter and the flux cutter at the turbine. The final test was con- 
ducted on 23 November 1966, when testing was halted due to time con- 
siderations. Under the present test rig set-up, the flux cutter and 
pickup is located inside the hood. It is felt that it would be advan- 
tageous to move the unit outside the hood where it would be readily 
accessible for adjustment. 

Data Reduction and Results 

The turbine flow rate was determined as described in the previous 
section. Total inlet temperature was determined from a probe installed 
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in the 6 inch pipe of the stator assembly. Inlet total pressure was 
taken from the average of the six Kiel probes (Figure 15) placed 
radially at the stator inlet. The hood pressure taken from the average 
of the eight hood static pressure taps was used for the turbine discharge 
static pressure. The turbine speed and dynamometer moment complete the 
quantities required to determine turbine performance. Data was reduced 
using the method of Vavra by a computer program. The method of analy- 
sis is being published in a Thesis by Lt. M. W. Wallace, U.S. Naval 
Postgraduate School, Monterey, California. [13] 

The following referred values were obtained from the performance 
analysis ; 





for 9" 


= cp/c.v = 1.4: 
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= Tto / 51 8.4 
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M ref 


= Mr = referred moment (Ft-lb) 


(52) 


N rtf 


= f\| r = referred turbine speed (RPM) 


(53) 


Wvtf 


=Wl/© = referred flow rate (lbm/sec) 
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(54) 


HPvcf 


= HP = referred horsepower (HP) 

VoJ 


(55) 


Ki s 


= = head coefficient 

UV2 



The turbine efficiency is a total to static efficiency, defined by 



(56) y\ * HP 

H Pth 

6 Vo 

where HP^ is the theoretical horsepower based on an isentropic expan- 
sion from the inlet total temperature to the static pressure in the 
hood, or, 



38 



(57) 



H Pi- Vi 
£ Vo 



-m) i' -(&) 



izi 

* 



703987 

(9-\ ) 



where 

P 2 = discharge pressure (lb/ft 2 ) 

Pro = inlet total pressure (lb/ft 2 ) 

Determination of turbine Reynolds numbers is rather arbitrary as 
to choice of characteristic lengths and velocities. After a perusal 
of the literature on the subject of Reynolds number effects in turbo- 
machines, two sets of quantities were chosen, one based on rotor tip 
speed, and one based on the flow rate and mean radius. They are de- 
fined as, 

(58) RE * W 

cAA Y 



and 

(59) RT * 2 TY RPM p R 

710 ^ 

where 

W = weight rate of flow (slugs/sec) 

2 

*aa = viscosity at turbine inlet total conditions (lb. sec/ft ) 
p = density at inlet total conditions (slugs/ft 2 ) 

R = rotor radius (inches) 

Yv* = mean rotor radius (ft) 

The referred values of horsepower, moment, and flow rates are shown 

in Figures 20 through 22. Turbine efficiency plotted as a function of 

the isentropic head coefficient is shown in Figures 23 and 24. Although 

great care was used in attempting to keep the pressure ratio at the set 

i 

value, the wide number of system variables caused minor deviations from 
the ratio. Table IV gives the actual pressure ratio at each data point 
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and the maximum deviation. 



4. Discussion of Results 

General 

The exhauster system functioned well, maintaining a hood pressure 
of about one half of the atmospheric pressure. At a pressure ratio of 
1.65 with a turbine flow rate of 3.88 lbm/sec a bothersome low fre- 
quency instability was noted in the hood pressure which was reflected 
in turbine operation. This fluctuation resulted in a low frequency 
flow rate fluctuation of as much as two inches of water in differential 
pressure across the flow nozzle and an RPM variation on the order of 
150. All of the turbine pressure readings were made on one manometer 
bank by using a Polaroid camera to photograph the board, and since the 
flow rate differential pressure, turbine RPM and dynamometer reading 
were taken as quickly as possible, the fluctuation did not seem to in- 
duce undue scatter in the data. 

The torque capsules are calibrated prior to each run by applying 
known weights to a lever arm. For dynamometer calibration, the torque 
capsule must be removed from the unit and placed in a special stand. 

The present calibration technique involves, jiggling or tapping the cap- 
sule as the weights are varied through several cycles until the read- 
out produces consistent values. For the last run, the stator axial 
force capsule could not be calibrated at all, but since data reduction 
by the method used does not require this value, no attempt was made to 
correct this difficulty. It is felt that the present method of cali- 
brating the torque capsules mechanically is unsatisfactory both from 
the standpoint of accuracy and time. Calibration of the torque cap- 
sules required as much as three hours, severely limiting the amount of 
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experimentation that could be accomplished during the normal working day. 
Provided that the capsules are not disturbed, electrical calibration 
would be both rapid and accurate. The present axial force measurements 
when used in a momentum analysis do not produce results consistent with 
those obtained from continuity, and the difficulty almost certainly lies 
with the axial force measurements. 

The turbine is designed for an isentropic head coefficient of 2.609 
at a pressure ratio of 1.5. Neither increasing the pressure ratio nor 
increasing the isentropic head coefficient did result in an increase in 
efficiency (Figures 23 and 24). At pressure ratios of 1.3 and 1.65 a 
dip in efficiency was noted with increasing head coefficients that is 
not present at the other pressure ratios. The reason for this dip is 
unknown, and it does not match the theoretical data given by Vavra. [in 
It is to be noted that this particular turbine is capable of operating 
under off-design conditions with very little reduction in efficiency. 

Although a great deal of research has been devoted to investigating 
Reynolds number effects in turbomachines, the results are still incon- 
clusive. Figure 25 shows the Reynolds number based on flow rate plotted 
versus turbine efficiency for three runs. Run 50 was conducted by the 
author using the full hood and exhauster installation. Runs 40 and 45 
were conducted by Professor M.H. Vavra with the turbine discharging to 
the atmosphere at pressure ratios of 1.3, 1.4, and 1.5. Run 45 used 
the same axial clearance of 0.41 inches and tip clearance of 0.033 
inches that were used on Run 50, whereas for Run 40, at an axial clear- 
ance of 0.410 inches the tip clearance was set to 0.015 inches. For the 
plotting of the data, the efficiency at each pressure ratio was averaged 
and plotted as a single point since the Reynolds number variation at a 
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given pressure ratio is insignificant. Reynolds number based on rotor 
tip speed is shown plotted versus turbine efficiency in Figure 26. 

NASA investigators on a number of turbines using the flow rate 
definition of Reynolds number found that turbines operating with Rey- 
nolds numbers above 2 x 10^ showed no variation in loss parameters with 
Reynolds numbers. f8] They also found that different turbines operating 
in the same range of Reynolds numbers had different loss parameters. 

This seems to point out that either the losses in turbines operated above 
a Reynolds number of 2 x 10~* are not a function of Reynolds number or 
that the machine characteristics used to determine the Reynolds number 
are not correct. 

Examination of Figure 26 which uses the Reynolds number based on 
rotor tip speed fails to show any trend, and again it appears that the 
efficiency is not a function of the Reynolds number within the range of 
investigation. 

Losses in turbomachines are due to viscous effects and since Rey- 
nolds number is a measure of these effects, it would seem only natural 
that some correlation should exist between Reynolds number and turbine 
efficiency. The primary difficulty lies in the elusiveness of the 
actual loss mechanisms associated with the flow in turbomachinery, and 
therefore, the inability of the investigator to be capable of relating 
the loss parameter to Reynolds number. A more detailed study would in- 
volve various definitions of characteristics lengths and velocities over 
a wide range of Reynolds numbers. This is a difficult task experimentally 
when dealing with a compressible fluid due to the wide range of tempera- 
tures and pressures required to obtain a significant variation in Rey- 
nolds number. 



42 



5. Conclusions and Recommendations 

The ARES MOD II turbine appears to reach its maximum efficiency at 
the design pressure ratio and the design value of the isentropic head 
coefficient. However, it is capable of operation over a wide range of 
pressure ratios and isentropic head coefficients with only a slight re- 
duction in efficiency. For the operating ranges tested, turbine per- 
formance does not seem to be related to Reynolds number. 

The weakest point of the present turbine test rig is the torque 
capsule calibration. A suitable means of electrical calibration should 
provide both speed and accuracy in calibration. The axial force 
measurements, both on the stator and coverplate, seem to be in error, 
due to some unknown cause. Therefore, it is recommended that the pre- 
sent instrumentation be improved, possibly by installing some type of 
hydraulic measuring device. 
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Table I 
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Long Mixing Tube 

Short Mixing Tube 

Subsonic Nozzle Operating Region 
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Table II 



Centerline Velocities, ft/sec, Pin/Pamb = 2.22 



St at ion 




Theoretical 


Experimental 


1 




1400 


1100+ 


2 


-* 


1270 


1260 


3 




920 


960 


4 




7 35 


760 


6 




570 


650 ++ 



+ Station 1 is the station where the centerline total pressure 
is much lower than the off-center total pressure, and the experi- 
mental value is in question. 

++ Station 6 is located sufficiently far downstream of the 
point where the mixing process is complete and the flow is that of 
turbulent flow in a pipe and is not expected to cbnform to the 
power law assumed for the centerline velocity- 
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Table III 



<J> and <J> cr itical as a Funct i° n of Secondary Mass Flow Rate 



Process from Secondary Flow Inlet to Annulus 



$> 


(^critical 


W(lbm/sec) 


Y\ 


.2614 


.2667 


3.43 


1.056 


.2755 


.2759 


3.62 


1.060 


.2876 


.2891 


3.785 


1.066 


.2986 


.3119 


3.93 


1.0770 


.3046 


. 3348 


4.00 


1.089 





Process from Hood 


to Annulus 




$ 


4>critical 


W(lbm/sec) 


Y) 


.3509 


.5950 


3.43 


1.295 


.3952 


.5874 


3.62 


1.287 


.4475 


.6015 


3.78 


1.302 


.5150 


.60 34 


3.93 


1.304 


.5616 


.6125 


4.00 


1 . 3141 
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TRAVERSING PROBES 
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TOTAL PRESSURE 
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FIG. 16 

TURBINE STATOk BLADE PROFILE 
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TUB BINE ROTOR AN'D STATOR DETAILS 
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